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a. 

Abstract 

The Klein-Gordon and Dirac equations in a semi-infinite lab (x > 0), 
in the background metric ds 2 = u 2 (x)(— dt 2 + dx 2 ) + dy 2 + dz 2 , are 
investigated. The resulting equations are studied for the special case 
u(x) = 1 + gx. It is shown that in the case of zero transverse-momentum, 
' the square of the energy eigenvalues of the spin- 1/2 particles are less 

i than the squares of the corresponding eigenvalues of spin-0 particles with 

f^**) . same masses, by an amount of mghc. Finally, for nonzero transverse- 

momentum, the energy eigenvalues corresponding to large quantum num- 
fS| , bers are obtained, and the results for spin-0 and spin- 1/2 particles are 

■ compared to each other. 

7' 

u 

bX). 



1 Introduction 



The behaviour of bosons and fermions in a gravitational field, has been of inter- 
est for many years, from the simplest case of a nonrelativistic quantum particle 
. in the presence of constant gravity, £Q for example, to more complicated cases 

5— i ' of the relativistic spin- 1/2 particles in a curved space-time with torsion, 0C3E] 

for example. Several experiments have been performed to test theoretical pre- 
dictions, among which are the experiments of Colella et al. |o], which detected 
gravitational effects by neutron interferometry, and the recent experiment of 
Nesvizhevsky et al. [H] , in which the quantum energy levels of neutrons in the 
Earth's gravity were detected. 

Chandrasekhar considered the Dirac equation in a Kerr-geometry back- 
ground, and separated the Dirac equation into radial and angular parts, 0|S]. 
In P|, the angular part was solved, and in JUj some semi-analytical results 
for the radial part were obtained. Similar calculations were performed for the 
Kerr-Newman geometry and around dyon black holes in refs. and ^2], re- 
spectively. 
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In this article, we investigate rclativistic spin-0 and spin-1/2 particles in a 
background metric ds 2 = u 2 (x)(— dt 2 + dx 2 ) + dy 2 + dz 2 , where the particles 
exist in a semi- infinite laboratory (x > 0). The wall x = 0, which prevents 
particles from penetrating to the region x < 0, corresponds to a boundary 
condition. For spin-0 particles, this is simply the vanishing of wave function 
on the wall. For the spin-1/2 particles, it is less trivial and will be discussed 
in the article. We consider the Hamiltonian-eigenvalue problems corresponding 
to a general function u(x), and obtain the differential equations and boundary 
conditions corresponding to the spin-0 and spin-1/2 particles. The special case 
u(x) = 1 + gx, is investigated in more detail. An exact relation between the 
square of the energy eigenvalues of spin-0 and spin-1/2 particles, whit same 
masses and no transverse momenta, is obtained, namely E^ = E^ G — mghc. 
Finally, the energy eigenvalues for large energies are obtained. 

2 Review of the non-relativistic problem 

The potential energy of a non-relativistic particle in a constant gravitational 
field is V^ rav (a;) = mgx, where g is the acceleration of gravity, the direction of 
which is along the x axis, towards the negative values of x. 

The Schrodinger equation is written as Htp^ch = *^<9t?/>s c h, where H = 
P 2 /(2m) + Vgrav, subject to the boundary conditions that ipsch does not di- 
verge at x — > ±oo. 

Writing ?/>sch = exp[(—iEt + ip2y + ip3z)/H]F(x), it is easily seen that F(x) 
must satisfy F"(x) = (L~ 3 x + L^ 2 Xs c h)F(x), where a prime means differen- 
tiation with respect to the argument, L _2 As c h := (p| +P§ — 2mE)/(h 2 ), and 
L := (2m 2 g/h 2 )~ 1 ^ 3 . This equation has two linearly independent solutions: 
the Airy functions Ai(L _1 a; + As c h) and Bi(Z a: + Asch) (see for example p. 
569 of [T3*|). Bi violates the boundary condition at +oo. So the solutions are 
Ai(Z _1 a: + As c h). These functions don't tend to zero at x — ► — oo, and it is 
expected, since the potential mgx is not bounded from below. But we note that 
a lab usually has walls. Consider a semi-infinite lab, for which 



In such a lab, the boundary condition at — oo is replaced with lim x _ > o+ V-'Sch = 0. 
So the solutions are Ai(L -1 a; + As c h), where As c h must be one of the zeros of Ai. 

The first four zeros of Ai are approximately Ai = —2.3381, A2 = —4.0879, 
A3 = —5.5206, and A 4 = —6.7867. If the transverse momenta of the particle 
vanish, the energy levels are these numbers multiplied by — h 2 ^ 3 (2m) 1 ^ 3 (g/2) 2 ^ 3 , 
the value of which for a neutron, in a lab on the Earth, where g = 9.8 ms -2 , 
is 0.59 peV. Therefore, the first 4 energy levels of a neutron are E\ — 1.4 peV, 
E 2 = 2.5 peV, £3 = 3.3 peV, and E 4 — 4.0 peV. This result has been recently 
verified experimentally 6 . 

3 A relativistic quantum particle 

According to general relativity, gravity is represented by a pseudo-Riemannian 
metric ds 2 = g^dx^dx" . We use the signature ( — h ++) for the metric. The 
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form of the metric depends on both the gravitational field and the coordinate 
system used to describe the field. 

In a spacetime with metric g^ Vl the Klein-Gordon equation, the equation for 
a spinless massive particle, is 



^^^ 99 d^'™ J^kg = 0, 

where g := det We have used a system of units in which the numerical 

values of the velocity of light (c) and the Planck constant (divided by 2ir) are 
both unity. 

To write the Dirac equation in a curved spacetime, or a flat spacetime but 
in a curvilinear coordinate system, one can employ the Equivalence Principle - 
see for example pp. 365-373 of |14j . but note that our convention for the Dirac 
Lagrangian, and therefore the Dirac equation, is different from that of Weinberg: 
we use '—to' in the following equation, instead of a '+to' in Weinberg's. It reads 
as follows: 

7° {da + r o ) Vd - mipD = 0, 
where r a s are spin connections, obtained from the dual tetrad e a , through 

de Q + r a b A e b =0 

T a 6 :=r a cfc e c 

T a :=-S bc T c a b 
Sbc ■=- 2 [76;7c] ■ 



We consider a gravitational field which is represented by the metric 



ds 2 = u 2 {x) (-dt 2 + dx 2 ) 
One can write the above metric, also like 



dy 2 



(2) 



ds 2 



-U 2 {X)dt 2 + dX 2 + dy 2 + dz 2 



where (dX)/(dx) = U(X) = u(x). In a small region of space (compared to 
the length in which the gravitaional field changes significantly) one can use 
u(x) = 1 + gx, U(X) = 1 + gX, and X = x. We will investigate the case of a 
general u(x), but then limit ourselves to the special case u(x) = 1 + gx. 
For the metric J5J, the Klein-Gordon equation reads 



dt 2 



dx 2 



dy 2 



dz 2 



IpKG = 0. 



To write the Dirac equation, we need the spin connections. The only non- 
vanishing spin connections for metric © are r°i = T 1 ^ = (u'/u 2 )e°. Thus, 
T 1 o° = u'/u 2 , and r = 7o7iu'/(2w 2 ). Therefore, the Dirac equation 



o 1 



i o 

for this metric reads as \^/ a d a 



2 7i — 



7 dt + 7 dx 



+ o7 — +7 7T+7% 
2 u A ay oz 



ipD = 0, which means 
d „ d 



i'B = o. 



(3) 
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4 Boundary condition at the infinite barrier 



In the non-relativistic equation, the infinite potential barrier which prevents 
particles from penetrating the region x < 0, leads to the boundary condition 
Um a ._ > o+ ^Asch = 0. This boundary condition emerges from the fact that the 
Schrodinger equation is second order in x, from which it follows that V'Sch must 
be continuous at x = 0. 

The Klein-Gordon equation is also second order in x, so the same boundary 
condition 

^kg(O) = (4) 

emerges. But the Dirac equation is of first order, so the four-spinor ?/>d can 
be discontinuous at x = 0, if the potential goes to infinity there. To find the 
proper boundary condition, one must first find a proper way of confining a Dirac 
particle to the region x > 0. The first guess is to add a step potential to the 
Hamiltonian. But this is the same as adding an electrostatic potential, which 
is the time component of a four-vector, having opposite effects on particles and 
anti-particles. Such a potential will not result in a wave function decaying as 
x — > — oo. A better way is to add a term — Vipu ipu to the Lagrangian, where 
V is a function of t, x, y, and z - something like a Higgs term. As a result 
of adding this 'scalar' potential to the Lagrangian, the mass m of the Dirac 
praticle is replaced by m + V. So the Dirac equation reads 

Now let's consider the step potential 

V( X ) = { V > X< ° (5) 
w [0, x > w 

where V is a positive constant. If the energy and the transverse momenta 
are finite, and the constant V is large, one can neglect the parts containing 
derivatives with respect to t, y, and z, arriving at 

( C l 1 d x -m-V{x))^{x) = 0. 

The solution to this equation, which does not diverge as x — > — oo, is V'd oc 
exp(qx) for x < 0, where q is a positive constant. As this equation is of first 
order, the solution must be continuous and we must have 

( 7 1 q-m-V r )^D(0 + ) = q 2 = {m + V) 2 =*> q = m + V (q>0). 

This means that 

( 7 1 -1)^d(0)=0, (6) 

which must be true, also for V — ► oo. Therefore, instead of Q, © is the 
boundary condition for a Dirac particle at the floor of the lab. 



5 The Klein-Gordon equation 

Since u in metric J2J does not depend on t, y, and z, one can seek a solution 
whose functional form is as ipKG{t,x,y, z) = exp(— iEt + ipiy + ip3z)ipKG{x)- 
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We arrive at 

H 2 . _ 1 

rhs.G(x)=0, x>0, (7) 



E 2 + ^-(p 2 + m 2 )u 2 (x) 



subject to the boundary condition V'kg(O) = 0, where p := y/ {P2) 2 + fe) 2 - 
Defining 

e := \/p 2 + to 2 , £ := ex, (8) 

we get 



d \ / d \ du 



6 The Dirac equation 

The functions in the equation J2J, do not depend on t, y, and z; therefore 
energy E and transverse momenta P2 and p^ are constants. By a suitable choice 
of coordinates (and without loss of generality) one can set p% = 0, and P2 = p. 
So we seek the solution as ipi)(t, x, y, z) = exp(—iEt + ipy)ijjjy(x) . Inserting this 
ansatz in Q , the resulting equation in terms of the a and j3 matrices becomes 

E i d i u' \ 
—/3 + -f3a 1 — + -j3a 1 ^-pf3a 2 -m)il> o = 0. (10) 
\u u dx 2 u z J 

Defining 

ijj := vV/'d, 

one arrives at 

E + ia%- upa2 — mu(3 ) ib = 0. (11) 

dx / 

In terms of the so called long and short spinors 

0:=i(l+j9)^, X:=\{l~P)h (12) 

which are eigenspinors of j3 with eigenvalues +1 and —1, respectively, (111(1 be- 
comes 

(E - mu) 4> + (^ ai ~^~ ~ upa^J X = 0, 

iai- upa2 ) 4> + (E + mu) x = 0. 

dx J 

Defining further \ := a\x, which satisfies (3\ — X> the above equations read as 
(E - mu) 4> + ~ upa 2 ui J X = 0, (13) 



i- upaia2 W> + (i? + mu) \ = 0. (14) 

da; / 

The matrix S defined through 

S := —ij3a\a2, (15) 
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commutes with a%, ct2, and /3; and its eigenvalues are ±1. So one can take ipu, 
and hence tjj, </>, \, and x, eigenspinors of S with the same eigenvalue in terms 
of which equations (|13|l and l|14|) can be written as 



(E - mu) i + qpu) 
- ?pu) (-B + mu) 



or 



X 



muoj, + spu<J2 — 



X. 

_d_ 

d.r 



0, 



where crs are the usual Pauli matrices. Defining 



e :=\Jp 2 + m 2 . 

:= arctan — , 
m 

we have 

ttkj 3 + <rpo2 = £ (03 cos 9 + 02 sin t 
Defining 



7T 7T 

2 2 



X 



one gets 



exp 



E 



X 



eexp -<Jip a^exp — —<j\ 



<p COS 2 ~ l X sm 2 

— z^sin I + xcos I 



EUtT3 — l(T\ — 

ax , 



) 






X. 



which can be written as 



E — eu 



dx 



i— E + eu 
ax 



= 0. 



Introducing <fi :— ± i%, this equation is transformed to 



E — eu 

dx 



da; 



£1* 



(16) 
(17) 



(18) 



(19) 
(20) 



(21) 



which, in terms of the variable £ = ex, leads to the following second order 
differential equation for <fi~ 



can be obtained from 



+ U)(f) 



- by 

e) E \d£ 



is 37 + u )<tr - 



(22) 



(23) 



Now turn to the boundary condition. In terms of <j> and \, eq. © is written 



as 



iX) U=o = 0. 



(24) 
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Using H19|l. we have 



iX = <t> cos - + 0+ sin-, 



from which and l)23|l . the boundary condition for tfi follows. 

r (o) cos °- + | (± + u^j r (0) sin °- = 0. (25) 
What remains, is to solve equation (|2~2")) subject to the boundary condition (J^SJ. 

7 The special case = 1 + gx 

7.1 The eigenfunctions 

For the special case u(x) = 1 + gx, we introduce 

X: =# + V? (26) 

and the Klein-Gordon equation becomes 



/kgPO = A| g /kgW, *>t/-, (27) 



where 

\2 ^KG 

A KG •— ) 

and 

/kgPO := V'kg (•£ ) • 
The boundary condition reads (see equation (J2J) 

/K G (^/f)=0. (28) 

For the Dirac equation, in this special case, again in terms of the variable 
X, equation i122[) becomes 



where 

E? 



\2 -^D 

An .— . 



and 

/ D (X) :=f(,). 
The boundary condition l|25|) . is now 

/d (vl) c °4 + £ (dl + *) /d (V!) sin I - °- (3o) 
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Comparing l|27|) and (|29[) . and defining A 2 = — 1 for the Klein-Gordon 
equation and A 2 = Aq for the Dirac equation, we see that both of them are of 
the same form: 

^ + x ){^ + x )^ = x ^ x >- < 31 > 

but subject to different boundary conditions l|28|) and H30|) . 

Note that this final equation is the same as the Schrodinger equation for a 
one-dimensional simple harmonic oscillator. The difference is that for the simple 
harmonic oscillator, — oo < X < +00; while in our case y/s/g < X < +00, and 

that we have the boundary condition l|28(l or (|30|l for X — > y / e/g + . 
To solve (|31|1 . we define h(X) as 

f(X) = h(X)cxp(-±X 2 ), (32) 



and obtain 

I 9X - 

dX J dX 

Now we seek a power-series solution for h(X): 

00 

h(x) = Y. a -- XT 


Putting this in one obtains 



«" = (33) 



(2fc)!r(~f) 



fl 2fe = ^~ / — — a o, (34) 



4 fc r (k + i - ^ 

fl 2fe+i = 77T; TTTTTTi PT ai ' ( 35 ) 



(2A + l)!r(|- 4 , 
and from that the following two soultions. 

°° (2X) 2fc r (fc- A 

r(-f)^ P*) 

=o (2X) 2fe+1 rffc + i- ^ 

We have to find a linear combination of these two functions, which remains finite 
as X — > 00. To do so, we define the function 

00 (2X) 2k+a T + f - x) 
Sa{X) ^ (2kTa)\ ' 

It is seen that ho and h\ are proportional to So and Si, respectively. Using a 
steepest-descent analysis to obtain the large- X behavior of S a , it is seen that 
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this behavior is in fact independent of a. This comes from the fact that if one 
replaces the summation over k with an integration, and change the variable 
k + a into k, then the dependence of S on a comes solely from the lower bound 
of the integration region. But this bound is unimportant, since the major part 
of the sum comes from large fc's. In fact, one can perform the steepest descent 
analysis and find 

S a (X) ~ X- 1 -^ exp(X 2 + l + ^). 
So, ho(X) + hi(X) remains finite as X — > oo, provided 

r(-f)- : ' 

So the unique (up to normalization) normalizable function which solves ()31|) , 

is 

h(X):=bJ2 f L - (38) 



7.2 Comparing the energy eigenvalues 

To compare the energies of three different Hamiltonians, namely the Schrodinger 
equation, the Klein-Gordon equation and the Dirac equation J2HJ), we 

reintroduce the physical constants previously chosen to be equal to 1. One 
obtains 

#KG 

E D 

To compare Ad and AkGj we first note that if p — 0, then 9 = (see equation 
(|18|l ). Therefore, the two boundary conditions H28fl and (|30|l become the same, 
and eigenvalue A 2 in (|31|l . becomes the same for the Klein-Gordon and the Dirac 
case. But A 2 is equal to X^ G — 1, and Aq. So, 

a d = a kg ^ !j for P = °! 



2m 



^KG 



Sch 



hg 



2/3 



= A D 



mc 3 



1/2 



mc 1 



mc , 

2 \ 1/4 



1/2 



mc" 1 + 



1/4 



or 

e d = e kg -mghc, for p = 0. (39) 

This is an exact result which determines the effect of spin on the gravitational 
interaction of relativistic particles. A fermion and a boson with same masses, 
have different quantum energies when fall vertically in a constant gravitational 
field. 

For p 7^ 0, the relation between the eigenvalues is more complicated, because 
although the differential equations are the same, the boundary conditions are 
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different. In this case, it is possible to find approximate relations for Akg and 
Ad, for large values of them, and compare them in this region. 

A WKB analysis of the wave function h(X), performed in appendix A, shows 
that 

A 2 



h(X) 



COS 7T XX 

4 



(40) 



Inserting 1(32(1 in ((30() . the boundary condition on h(X) for Dirac particles reads 

1 dh . 



'e/g 



= 0, 



which, using (|40|l . leads to 



7rA 



D XnX 



X=Je/g 



or 



7rA D !~e 6 f 1 



(41) 



For the Klien-Gordon equation, from l|28[) and 1)32(1 . it is seen that h{^e/ g) 
0. Using 1)40(1 . and remembering A 2 = Aj^Q — 1, we find 



7T, 



(42) 



in which 0(1/Akg) terms have been ignored. Comparing ((41(1 and 1(42(1 . results 
in 

26» 



— ^KG - H > 
7T 



or 



e d = e kg ~ Vp 2 + m2 c 2 ^ 5 ( 1 ~ ~ arctan 

7T 



m 



(43) 



A The asymptotic behaviour of the function h(X) 

To obtain the asymptotic behaviour of h(X), ea. l(38|) . for A 3> 1 and finite X, 
we first write ((38|l as 

h(X) = bY , \ 4 -l-bY - V ' ^. (44) 

fro ( 2fc ) ! fro ^ + 

Using r(p)r(l — p) = ir/ (sinp7r), we have 



4; dn^fc-^r^-fe + l) sin(7r^)r(^-fc + l)" 

(45) 

As X is finite, large fc's have negligible contributions in the power series for 
h(X). So we can use the Stirling's formula T(x + 1) ~ ^/2nx x+ ^ 1 / 2 '' l e~ x , to 
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obtain 



r — - k + 1 =V2n[ — 



2X 



k 



4 



A 2 /4 

-A 2 /4 



— fcH 



.fe-(A 2 /4) 



Therefore (|45|) leads to 

rifc-- 



(-1)*tt ^A 2 ^- (1/2 V4e^ A2/4 

A 2 



%/2^sin(7r^) V 4 



A similar argument shows that 



(-1)*7T 



27TC0S(7T^) V 4 



4e\ 
A 2 "] 



Inserting JJBJl and (|4*7jl in (|4*Ijl . one obtains 



h(X) 



,/A 2 ^ 1/2 



2 V 4 



A 2 /4 



(46) 



(47) 



1 



-l) fc (AX) 



2/r 



1 °° 



(~l) fc (AA) 
(2* + l)! 



2fe+l 



2bV2iT 



4e 



A 2 /4 



A 2 

cos | 7r AX 

4 



Asin(Trf) \X\ 
This is nothing, but the leading-order result of a WKB analysis. 



(48) 



11 



References 

[I] L. D. Landau and E. M. Lifshitz; Quantum Mechanics - Non-relativistic 
Theory (Pergamon Press, 1977). 

[2] D. V. Ahluwalia & C. Burgard; Gen. Rcl. Garv. 28 (1996) 1161. 

[3] C. Y. Cardall & G. M. Fuller; Phys. Rev. D55 (1997) 7960. 

[4] M. Alimohammadi & A. Shariati; Int. J. Mod. Phys. A15 (2000) 4099. 

[5] R. Colella, A. W. Overhauser, & S. A. Werner; Phys. Rev. Lett. 34 (1975) 
1472. 

[6] V. V. Ncsvizhcvsky et al; Nature 415 (2002) 297. 

[7] S. Chandrasckhar; Proc. R. Soc. Lond. A349 (1976) 571. 

[8] S. Chandrasckhar; The Mathematical Theory of Black Holes (Oxford Uni- 
versity Press, 1992). 

[9] S. K. Chakrabarti, Proc. R. Soc. Lond. A391 (1984) 27. 

[10] B. Mukhopadhyay & S. K. Chakrabarti; Class. Quant. Grav. 16 (1999) 
3165. 

[II] D. N. Page; Phys. Rev. D14 (1976) 1509. 
[12] I. Scmiz; Phys. Rev. D46 (1992) 5414. 

[13] C. M. Bender & S. A. Orszag, Advanced Mathematical Methods for Scien- 
tists and Engineers (McGraw-Hill, 1978). 

[14] S. Weinberg; Gravitation and Cosmology; (John Wiley, 1972). 



12 



